Localization and Mobility Edge in One-Dimensional Potentials 

with Correlated Disorder 
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We show that a mobility edge exists in ID random potentials provided specific long-range cor- 
relations. Our approach is based on the relation between binary correlator of a site potential and 
the localization length. We give the algorithm to construct numerically potentials with mobility 
edge at any given energy inside allowed zone. Another natural way to generate such potentials is to 
use chaotic trajectories of non-linear maps. Our numerical calculations for few particular potentials 
demonstrate the presence of mobility edges in ID geometry. 
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It is commonly believed that there is no mobility edge 
in ID models with random-like potentials. This is based 
on the fact that for random potentials all eigenstates 
are exponentially localized, no matter how weak the ran- 
domness is 0. On the other hand, for potentials with 
"correlated disorder" the localization length diverges for 
specific values of energy (see, e.g. ^,||)- The well-studied 
model of this kind is the so-called "random dimer" [Q for 
which the potential has peculiar short-range correlations. 
Though there is no mobility edge for such potentials, this 
example shows a highly non-trivial role of correlations. In 
this Letter we study the relation between correlations in 
the site-potential of ID tight-binding model and localiza- 
tion properties of eigenstates, and give examples of the 
potentials with mobility edges inside the energy band. 

The model under consideration is the discrete 
Shrodinger equation for stationary eigenstates i/} n (E) , 



where 



tp n +i + r/> n -i = (E+ e n ) ip n 



(1) 



where E is the eigenenergy and e n is the site-potential. 
To study the origin of delocalized states in long- 
correlated random potentials, we suggest a simple and 
clear approach based on the representation of the quan- 
tum model (|]J) in terms of classical two-dimensional 
Hamiltonian map, 



Pn+1 = Pn + (E - 2 - e„) x n , 

*^n+l — ' T" Pn+l 5 



(2) 



where p n +i = x n +i — x n and x n = ip n . This map de- 
scribes the behavior of a linear oscillator subjected to 
linear periodic delta-kicks with amplitude depending on 
e n . In such an approach, localized quantum states corre- 
spond to trajectories which are unbounded in the classi- 
cal phase space (p, x) when n — > oo. Contrary, extended 
states are represented by bounded trajectories. 

It is convenient to introduce the action-angle variables 
(r, 9) and represent Eqs.([l]) and (|J) in the followind form 
(see details in H), 



sin 

COS I 



'71+ 1 

'n+1 



D-^sm (0„ 
D- 1 (cos (0„ 



fi) - A n sin 6 n sin fi) , 
/i) + A n sin n cos ji) . 



(3) 



D n = ^ = V 1 + A n sin (20„) + Al sin 2 #„ 



A, 



-e„/sin/z, E = 2 cos /j . 



(4) 



We use the following definition of the inverse localization 
length (or Lyapunov exponent A) || , 



1 



N 



r 1 = A = lim — > ' In 

71=1 



'n+l 



(5) 



which coincides inside the energy band, \E\ < 2 with 
the standard definition ||, A = (In | ^ +1 1) (see details 
in Q). Here the brackets stand for the average over n. 
This Hamiltonian map approach turns out to be very ef- 
fective in the study of completely disordered potentials 
as well as potentials with correlated disorder ||[| ■ 

Below we consider a general case of any stationary site- 
potential e n under the condition \e n \ << 1. In this case 
one can expand the logarithm in Eq. (^|) and in the sec- 
ond order of perturbation theory gets, 



r 1 



(<£) (e„sin(20 n )) 



!sin 2 /i 



2 sin /i 



(6) 



In order to calculate the correlator (e n sin (20 n )) 
with quadratic accuracy, we use the approximate one- 
dimensional map for the phase n obtained from Eq. (^) , 



~>n — Dn-1 



sin fi 



(7) 



From Eq. (Q) the correlator (e„ sin (2#„)) can be ex- 
pressed through the precedent one, (e„ sin 2(0„_i)). Let 
us introduce the following notations, 



Q = ( p 2i6„- k \ _ i cxp(-2^) 2 

(e„e„_ fe ) = e^(fe) , (e 2 ) = e 2 . 



(8) 
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It can be shown that the normalized correlators a k are 
defined by the infinite set of linear equations, 



1 < k < oo. 



(9) 



These equations emerge after multiple application of the 
recursion relation (Q) to the correlator Q k . 

According to Eq. (||) the localization length is deter- 
mined by the only term ao, 



r 1 



-Re (e" 2v 



!sin 2 /z 4sin 2 /x 
which can be easily obtained from Eq. 



a = £(fc) exp [— 2i^i(fc — 1)] . 



(10) 



(11) 



fc=i 



As a result, we come to the final expression for the inverse 
localization length, 



i _ 4 



r 1 = 



i sin 2 /i 



¥?0) = l + 2]T£(fc) cos(2^fc) . (12) 



fc=i 



Here, the function ip (ji) is given by the Fourier series with 
the coefficients £(fc) which are the correlators of the site 
potential e n . Note that Eq. ( fl2"| ) fails close to the band 
center, E = 0, and to the band edges, 25 = ± 2 , where 
even the standard Anderson model exhibits peculiarities 
(see, for example, M and references therein). 

Eq. (|l^) establishes the relation between the localiza- 
tion length and correlations. It allows us to calculate 
the localization length if statistical properties of the se- 
quence e n are known. Note that only binary correlators 
enter in Eq. ([l2|). This property comes from the Born 
approximation and not from restrictions for statistics of 
e„ ||. In particular, we do not assume that the statistics 
is gaussian. 

Formally, the same result has been derived by differ- 
ent methods in ||[l(J. For continuous models, the high- 
energy asymptotics of the Lyapunov exponent is also de- 
termined by the correlation function of the site-potential 
|pL To the best of our knowledge, the relation ( |T2| ) has 
not been discussed yet concerning the problem of mobil- 
ity edge. The advantage of the above approach based on 
the Hamiltonian representation, is that it allows a simple 
generalization to the Kronig-Penney model with corre- 
lated disorder JllJ and to the models with off-diagonal 
disorder. The approach also clarifies the physical mean- 
ing of the summation in Eq. (|l2|) . Namely, the localiza- 
tion appears as a result of multiple scattering at different 
sites, see Eq.(||), and each scattering enters in the Born 
approximation, see Eq.(Q). 

To reveal the role of correlations, let us apply the gen- 
eral relation ( |T2| ) to some known models. First, we con- 
sider the potential 



y2 cos(27rcm 7 ) 



(13) 



with a irrational. The correlation function £(fc) for this 
sequence is given by the sum 



1 



N 



£(k) = lim — N ' cos(27T7a£;n 7 x ) 

71 = 1 



(14) 



For < 7 < 1 all the correlators £(fc) are independent 
on k and equal 1 , thus, giving ip (/x) = for < // < tt. 
Therefore, for this case the result coincides with the case 
of the constant potential e n = eo; in both cases the in- 
verse localization length vanishes. This is consistent with 
the perturbative approach developed in fl0| . 

For 7=1 Eq. ( l]) is reduced to the Harper equation 
with incommensurate potential. In this case the correla- 
tion function £(fc) oscillates, £(fc) = cos(27rafc), that also 
gives ip (/i) = 0. Thus, in the Harper model with a weak 
potential all states are extended [0,[L3|. For 7 > 1 the 
correlators £(fc) vanish and the localization length turns 
out to be the same as for the Anderson model (for energy 
close to the band center, see 

Let us consider now the random dimer model Q which 
is specified by the sequence e n having only two values e\ 
and €2, each of them appears in pairs (e\€\ or €262) and 
each pair emerges with probability 1/2 . Fully transpar- 
ent states are known to occur for E = e\ and E = ei , 
see, for example, || and references therein. Statistical 
properties of this model are characterized by the vari- 
ance e 2 = I (e 2 + e|) and two correlators, (e„e n _i) = 

,2 3+2A+3A 2 
fc 4(1 + A 2 ) 



and 



e 2(l + *i) wllC1 ' C A = 61 A 2 



By substituting these correlators into Eq. (12), we obtain 
the inverse localization length for the dimer, 



l-\E) 



ei(l-Ar 



E 2 



,(1 + A 2 ) 4-£ 2 



(15) 



The generalization for A^-mer when the values ei and 
62 appear in randomly distributed blocks of length N , 
can be readily done. For example, for the trimer we have: 



l-\E) 



egq-A) 2 (E*-iy 

12(1 + A 2 ) 4-E 2 



(16) 



Note that the above explicit expressions for dimer jl5| ) 
and trimer ( |l6|) are given for the whole energy range; to 
the best of our knowledge, this is a new result. 

For the N— mer constructed from the standard An- 
derson model by repeating N times each random value 
e n , Eq.(p^) gives the same result obtained in |l9]. In 
particular, this result shows that the correlations do not 
necessarily suppress the localization, they can make it 
even stronger than in the Anderson model [M. 

The localization length ([l2]) does not depend explic- 
itly on the site-potential but on the binary correlators 
£(fc)of the potential. Let us consider different types of 
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correlations. First, we take a sequence e„ with exponen- 
tial decay of correlations, = exp(— /3k) . This type 
of correlations is known to occur in completely chaotic 
Hamiltonian models without stable regions in the classi- 
cal phase space. By substituting £(fc) into Eq. (|l2|), one 
gets the following expression, 



r 



sinh j3 



8sin 2 /i cosh/3 — cos(2^i) 



(17) 



which establishes a link between the correlation radius 
(parameter /3 _1 ) of the site-potential, and the localiza- 
tion length of eigenstates. In the limiting case (3 — > oo, 
we get the result for uncorrelated sequence (Anderson 
model) and for (3 — > the result for the constant poten- 
tial is recovered. 

Concerning the power decay of correlations, one should 
note that such simple functions as £(fc) = 1/fcor 1/fc 2 can 
not serve as the correlation functions of any random se- 
quence. Indeed, these functions lead to negative values of 
the Lyapunov exponent and do not satisfy the condition 
required for any correlation function j^j. Specifically, 
the matrix = — k') must be positive-semidefinite, 
the property which is not obvious a priori for arbitrary 
positive function 

An important application of Eq. ( |l2| ) is related to the 
question about existence of the mobility edge in ID sys- 
tems. So far the mobility edge has been predicted for 
incommensurate potentials Eq. (|l3| ) with 7 < 1 and finite 
amplitude eo [f0| 19| , and for Kronig-Penney model with 
constant electric field pofl . 

The relation (|l^) allows us to construct site potentials 
with mobility edges. Rewritten in the form 



tt/2 

^(k) = — / if (fi) cos(2fc^) dfi 
n J 



(18) 



this relation gives a solution of the "inverse problem" , 
namely, it shows how to calculate the correlation func- 
tion £(fc) of the random potential if the normalized Lya- 
punov exponent Ao (E) = |</? (arccos(i?/2)) is known. In 
other words, for any dependence Ao there exist a set of 
sequences e„ with correlation function £(fc). 

To reconstruct the sequence via the correlation func- 
tion, we have used the following algorithm pjj, 



e„ = 4 ^ £(k)Z n+k , 



(19) 



k— — c 



where A is the normalization constant and Z s are ran- 
dom numbers from the interval [0,1]. As an example, 
we consider two sequences e„ which exhibit the mobil- 
ity edges at E = ±1. The first one was obtained from 
Eq.Q by substituting = affe sin(^) and A = 1. 
This corresponds to the step-function dependence of the 



Lyapunov exponent, Ao = for \E\ < 1 and Ao = 1 for 
1 < \E\ < 2. Circles in Fig. 1 show the dependence 
Ao(-E) calculated numerically for the reconstructed se- 
quence e„ using Eq. ([5]) . The step-function is reproduced 
quite well and oscillations at \E\ > 1 are mainly due to 
the finite length of the sequence, N = 10 6 . 




FIG. 1. Rescaled Lyapunov exponent Ao(E) for two 
site-potentials e n constructed by Eq. (^|) for 300 correlators 
£( fc ) = and n = l, 10 6 . 

Another sequence e„ was generated in order to get 
a smooth increase A (E) for \E\ > 1, see squares in 
Fig. 1. To do this, we used the correlators £(k) which 
are the Fourier components (|l8|) of the linear function 
Ao(-E') = 1.53(|£| - 1) for \E\ > 1. As one can see, ac- 
tual dependence Aq(E) in Fig. 1 for E > 1 clearly differs 
from the linear one. This results from the fact that in 
general case the algorithm (|l{|) reconstructs the sequence 
e„ with correlations £(fc) only asymptotically, for k 3> 1 
pl| . Nevertheless, it allows us to generate sequences with 
different values of critical exponents at the mobility edge. 
These examples demonstrate the existence of mobility 
edges for ID site potentials with long-range correlations. 

Another way to obtain explicitly the sequence e„ with 
a mobility edge is to use weakly chaotic non-linear maps. 
Let us take the well-known standard map |22fl , 



P n +i = P n + K sm2%X n (mod 1) 
X n+ i = X n + P n+ i (mod 1) 



(20) 



For K < 1 there are trajectories close to separatrices 
of non-linear resonance, with a very slow decay of cor- 
relations. For one of such trajectories (see the inset in 
Fig. 2) we have constructed the sequence e„ as follows, 
e n = 2eo sin(27rA^ n ), and calculated the Lyapunov expo- 
nent A (-E) = 8Asin 2 /i/e 2 , from Eq.(^). The dependence 
given in Fig. 2 for eo = 0.1 clearly indicates the presence 
of two mobility edges. It is important that the mobility 
edges are stable with respect to rather wide variations of 
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e . Independent computation of eigenenergies ofEq. © 
with the corresponding potential e n demonstrates that 
the energy region shown in Fig. 2, indeed, belongs to the 
energy spectrum. 




FIG. 2. Dependence Ao(E) for the model (g) with the site 
potential defined by Eq.(EOI). The length of the sequence 
e„ is N — I0 6 . The corresponding trajectory (Xn, Pn) for 
X = 0, Pn = 0.292 and K = 0.8 is shown in the inset. 



It is interesting to note that inside a narrow chaotic 
region like that shown in the inset of Fig. 2, the time 
dependence of trajectory (X n , P n ) is given by regular 
rotation around the resonance of period 3 with the fol- 
lowing chaotic motion in the vicinities of the crossed sep- 
aratrices. This type of motion can be compared to the 
so-called intermittency ||| which is well studied in ID 
maps. Therefore, one can expect that there is a direct 
link between the correlations in the site potential and 
the intermittency in the evolution of the corresponding 
dynamical model which generates this potential (see also 

6). 

In conclusion, we have shown the existence of the mo- 
bility edge in ID random potentials. The relation jl^ ) 
allows us to determine binary correlation function of the 
site potential which exhibits any desirable dependence 
of the localization length on energy. Specific problem we 
address in this Letter, is how to reconstruct the potential 
itself from the binary correlations. We suggest an algo- 
rithm to calculate it numerically, and present the data for 
few specific potentials, with two mobility edges inside the 
allowed energy zone. Potentials exhibiting mobility edges 
can be also constructed by making use of weakly chaotic 
trajectories of non-linear maps, this is demonstrated us- 
ing the standard map. Although analytical results are de- 
rived in quadratic approximation with random potential, 
our numerical data clearly show that the mobility edges 
do not disappear for relatively strong potentials. This 



fact opens a possibility for customer design of superlat- 
tices with anomalous electronic transport, in particular, 
devices with selective transparancy. 

Very recently, the so-called self-similar potentials with 
specific long-range correlations have been studied [p5[ , 
for which numerical data cleraly indicate the exictence 
of mobility edges close to the center of the enery band. 
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